%% Fig. 1H - R(T) for Cooldown.
figure(1);
hold off;
plot(RT.TSample, RT.L1 / RT.L1(1), 'linewidth', 2, 'color', [63,72,251]/255);
hold all;
plot(RT.TSample, RT.L2 / RT.L2(1), 'linewidth', 2, 'color', [249,61,96]/255);
plot(RT.TSample, RT.L3 / RT.L3(1), 'linewidth', 2, 'color', [104,170,251]/255);
xlabel('T (K)');
ylabel('R / R(290K)');
ylim([0, 1]);
xlim([50,280]);

%% Figure 1I - I-V in both directions.
figure(1); clf
hold off;
%8.7K data. Junction Area = 376 um^2; Voltage Preamp gain = 10 x. 
plot(IVf.K2400aA / 376 * 1e5, IVf.DMM1 / 10 * 1000 - 0.21025, 'linewidth', 2, 'color', [78,49,245] / 255);

hold all;
plot(IVr.K2400aA / 376 * 1e5, IVr.DMM1 / 10 * 1000 - 0.21025, '--', 'linewidth', 2, 'color', [78,49,245] / 255/2);
xlabel('Current Density (kA / cm^2)');
ylabel('Voltage (mV)');
xlim([-1.5, 1.5]);
ylim([-20, 20]);


%% Figure 1I Inset
figure(1);
hold off;
%8.7K data. Junction Area = 376 um^2; Voltage Preamp gain = 10 x. 
plot(IVf.K2400aA / 376 * 1e6, IVf.DMM1 / 10 * 1000 - 0.21025, 'linewidth', 2, 'color', [78,49,245] / 255);

xlabel('I / A (\muA / \mum^2)');
ylabel('V (mV)');
xlim([5, 13]);
ylim([0, 1.2]);
